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Abstract. Entropic aging consists in a progressive slowing down of the low- 
temperature dynamics of a glassy system due to the rarefaction of downwards 
directions on the energy landscape, as lower and lower energy levels are reached. 
A prototypical model exhibiting this scenario is the Barrat-Mezard model. We 
argue that in the zero-temperature limit, this model precisely corresponds to 
a dynamical realization of extreme value statistics, providing an interesting 
connection between the two fields. This mapping directly yields the long-time 
asymptotic shape of the dynamical energy distribution, which is then one of 
the standard extreme value distributions (Gumbel, WeibuU or Frechet), thus 
restricting the class of asymptotic energy distributions with respect to the original 
preasymptotic results. We also briefly discuss similarities and differences between 
the Barrat-Mezard model and undriven dissipative systems like granular gases. 



PACS numbers: 02.50.-r, 05.40.-a, 64.70.P- 

1. Introduction 

Glassy systems exhibit a slow relaxation towards their equilibrium state, giving rise to 
non-stationary and aging properties [T . Different physical mechanisms have been put 
forward to account for this intriguing phenomenon. One such mechanism is related 
to the presence of energy barriers that the system needs to overcome through thermal 
activation in order to explore its phase space. The system is then trapped for long 
times in energy wells, and successively explores deeper and deeper traps In such 

a scenario, thermal activation plays a crucial role and the dynamics becomes frozen 
in the low temperature limit. Another standard, though also schematic, mechanism 
leading to the aging phenomenon in glassy systems is the so-called 'entropic aging', 
where the dynamics is slowed down by the rarefaction of downwards directions. Energy 
barriers are then replaced by entropic barriers, and the temperature no longer plays a 
prominent role in this case. Explicit realizations of this scenario have been proposed 
in simplified models, like the Backgammon model [HIS], or the Barrat-Mezard model 
[6l [7] , the latter probably being the simplest model one could think of to illustrate the 
'entropic aging' phenomenon. 

In this short note, we establish a direct and simple link between the Barrat- 
Mezard model and extreme value statistics, namely the statistics of the maximum 
or minimum value in a set of random variables |S1 [5] • This link then allows for the 
direct application of the standard results of extreme value statistics to the Barrat- 
Mezard model, providing in a straightforward way some non-trivial information on 
the asymptotic energy distribution reached dynamically in the long-time limit. We 
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also discuss, in a more exploratory way, a possible qualitative analogy between the 
Barrat-Mezard model and the free relaxation of dissipative systems like granular gases. 

2. The Barrat-Mezard model 

As mentioned in the introduction, the Barrat-Mezard model [6| is a simple stochastic 
model exhibiting the aging phenomenon, due to the presence of entropic barriers. In a 
continuous energy formulation, valid in the limit of an infinite number of microstates, 
the model is defined as a continuous-time Markov process in which the system jumps 
from an energy E to another energy E' with a Glauber rate: 

where Tq is a microscopic frequency scale, p{E) is the density of microstates at energy 
E, and T is the heat bath temperature. The probability P{E, t) to be in a microstate 
of energy E at time t obeys the following master equation: 

^{E,t)^ p dE'[w{E\E')P{E\t)~w{E'\E)P{E,t)]. (2) 

The transition rate ([1]) satisfies detailed balance, so that P{E^ t) is expected to 
converge when t — oo to the equilibrium distribution 

P.^{E,t)^^p{E)e-^/^. (3) 

Depending on the shape of the density of states p{E), the equilibrium distribution 
Pcq{E,t) may not be normalizable below a given temperature meaning that the 

continuous energy description breaks down (the probability measure concentrates on 
the few lowest energy states). For instance, if p{E) is exponential, 

piE)^T-'e^/^^e{-E) (4) 

(with O the Heaviside function), the distribution Pcq{E,t) is non- normalizable for 
all temperature T < Tg. The slow relaxation (or aging) regime can however be 
described through a continuous energy description [7,. If p{E) decays faster than 
exponentially (e.g., if p{E) has a Gaussian tail), the equilibrium distribution is well- 
defined in a continuous energy framework for T > 0, and aging is interrupted beyond 
the equilibration time. 

An interesting property of this model is that, contrary to what happens for trap 
models involving thermal activation [31 13] , the dynamics is not completely frozen in 
the zero temperature limit. The average trapping time r needed to escape a given 
microstate at energy E is given for T = by 

-I /"OO t^E 

-=/ dE'w{E'\E)= dE'p{E'). (5) 

^ J — oo J — oo 

Microscopic configurations can then be labelled by r rather than by E, and the 
dynamics of the system is described by the probability p(r, t) to occupy a configuration 
characterized by t at time t. It has been shown [6j that at large time, the dynamical 
distribution p{t, t) takes a scaling form 

p{r,t) = j<p[^), (6) 
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with a scaling function given by 

1 



(7) 



This shape of the scahng function (O is generic and does not depend on the shape of 
the density of states p{E). However, the dependence on the functional form of p{E) 
is hidden in the definition of r given in Eq. ([5]). Indeed, expressing P{E,t) as 

(8) 



P{E,t)^p{T,t) 



dE 



one clearly sees from Eqs. ([5]) and (jB]) that P{E,t) depends on the detailed shape of 
p{E), suggesting that Eq. ^ only corresponds to a preasymptotic regime. 

Let us also mention that the scaling form ^ is not necessarily restricted to the 
zero-temperature limit. In the specific case of an exponential distribution p{E), as 
given in Eq. ([4]), a similar scaling form exists in the whole aging regime T < Tg, with 
a temperature-dependent scaling function (j) [7]. 



3. Mapping onto extreme value statistics 



It is interesting to observe that the above zero-temperature dynamics can be 
interpreted as a dynamical realization of extreme value statistics. To make the analogy 
more precise, let us first define the discrete time version of the Barrat-Mezard model. 
The dynamics is defined in the same way as above, apart from time discretization: at 
each time step, a new microstate of energy E' is picked up from the distribution p{E), 
and the move from E to E' is accepted with probability 

and rejected otherwise. Consistency between the discrete and continuous time versions 
of the model implies that the time step tq in the discretized version should be chosen 
as To = Fq"^, where Fq is defined in Eq. ([T|). 

It is then easy to see that the zero-temperature discrete-time Barrat-Mezard 
model exactly implements, in a dynamical way, an extremal process. Let us assume 
that at time i = 0, the initial energy Eq is drawn at random from the distribution 
p{E), which can be interpreted as an infinite temperature initial condition (namely, 
all microstates have the same weight). Then at time t = utq, the occupied value E{t) 
is the lowest value among the set of n + 1 random variables Ej, j = 0, . . . ^n, that have 
been drawn, and either accepted or rejected, at each time step: 

E{t)=rRm{Eo,Ei,...,E„). (10) 

Such a mapping from the Barrat-Mezard model to extreme value statistics is 
interesting for several reasons. First, it provides an example of a non-trivial dynamical 
process which maps exactly onto an extreme value problem. Second, standard results 
from the field of extreme value statistics can then straightforwardly be applied to 
the Barrat-Mezard model, giving a direct access to the long-time asymptotics of the 
distribution P{E,t). Indeed, it is known from extreme value statistics that the shape 
of P{E, t) asymptotically attains a universal form belonging to three possible classes, 
depending on the tail of the distribution p{E) [8l [9] . These universal distributions are 
obtained through an asymptotic rescaling, 

P{E,t)^^^<^>{^^^, f>To, (11) 
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with suitably chosen time-dependent parameters at and bt- If P^E) decays for 
E — !■ —oo faster than any power law, then the asymptotic distribution $(z) has a 
Gumbel (or Fisher-Tippett-Gumbel) form, namely 

$(z) = ^ e^-^\ (12) 

Note that we consider here the statistics of minimum values, so that the definition (fT2|l 
of the Gumbel distribution differs by a sign reversal from the more usual definition 
related to maxima. 

If p{E) instead decays as a power law p{E) ^ {a > 0) when E — oo, 

then ^{z) is the Frechet distribution 

$(z) = $f(z)^^eH-l"° (z<0). (13) 

Finally, in the case when E is bounded from below, and p{E) behaves as a power law 
close to the bound i?min, 

p{E) ^ {E ^ E^i,,r-\ -B^£;„,in iE>E^i^) (14) 

with J/ > 0, the asymptotic distribution is the WeibuU one 

^z) = ^„{z) = uz"-^ e-^" (z>0). (15) 

These three cases can be encompassed into a single distribution $7(2:), using a slightly 
different choice of the rescaling parameters at and bt }10| 

$7(2) = (1 - 7z)"^"^ exp -(l-7z)"^ , l-7z>0, (16) 

with a real parameter 7. By continuity, $7(2) should be interpreted for 7 = as 
exp[— exp(z)], which corresponds to the Gumbel distribution. For 7 > 0, Eq. ([T6|) 
corresponds to the Frechet case, with 7 = 1/a. For 7 < 0, the WeibuU distribution is 
obtained, with 7 = —l/iy. 

As a result, we have shown using a mapping onto extreme value statistics that 
the dynamical energy distribution P{E, t) of the Barrat-Mezard model asymptotically 
corresponds to one of the standard extreme value distributions (Gumbel, Weibull or 
Frechet), depending on the behaviour of the tail of the density of states p{E). This 
result significantly restricts the class of allowed asymptotic distributions with respect 
to the original results obtained by Barrat and Mezard. As mentioned at the end of 
Sec. [21 the distribution P{E, t) deduced from the scaling form ^ actually corresponds 
to a family of distributions indexed by the full density of states p{E), once the variable 
r is expressed as a function of E. In constrast, the present distributions obtained from 
extreme value statistics correspond to a one-parameter family of distributions, that is 
to a much more restricted class. This difference can be ascribed to the fact that the 
scaling regime (jB]) can be considered as a preasymptotic regime. 

The above results are valid, strictly speaking, in the case of discrete time 
dynamics. However, the type of dynamics (discrete or continuous time) should not 
influence the large-time asymptotic dynamics of the system, so that extreme value 
distributions are expected to hold also for continuous time dynamics. This is confirmed 
by the simple example of the exponential distribution (|4]), which falls into the Gumbel 
class. In this specific case, one has t{E) = e~^/"^», and the scaling form ^ with (j) 
given by ([7]) directly boils down to the Gumbel distribution. 
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4. Illustration on large systems 

In this section we briefly illustrate the above results on two different classes of 
systems: mean- field models like the Random Energy Model (REM) [TTl [1^, and 
finite-dimensional models with short-range correlations. 

In the REM, which can be considered as the simplest disordered system, all 
microscopic configurations have an energy E drawn at random from a Gaussian 
distribution p{E) with a variance proportional to N Given that the low- 

temperature properties of the model depend on the low-lying energy states, only the 
low energy tail matters. It has thus been proposed [T^] to generalize the Gaussian 
distribution to a distribution p{E) with a tail 

p(i?)^e-l^l'/^'", E^-oo (/3>1), (17) 

where iV ^ 1 is the underlying number of degrees of freedom. Algebraic prefactors 
can also be included, but they do not influence the asymptotic behavior [12]. The 
REM exhibits interesting equilibrium properties, like a freezing transition at a finite 
temperature Tc = /3-Hln2)(i-'^)/'' [T2^. 

In its original formulation, the REM is an equilibrium model, without any 
specified dynamics. A rather natural choice is to use a Glauber dynamics, connecting 
all microscopic configurations with the probability (|9]). It has been shown that for a 
finite number of states and a non-zero temperature, the dynamics of the REM becomes 
similar to that of the trap model beyond a size-dependent time scale, which diverges 
with system size TJ'. However, in the zero-temperature and infinite size limits, the 
dynamics by definition corresponds to the Barrat-Mezard model, with the density of 
states p{E) given in Eq. p7)) . From the results of Sec.|3l the corresponding long-time 
energy distribution P{E,t) is a Gumbel distribution, for all (3 > 1, since the density 
of states p{E) decays faster than any power law. 

As a second example, we consider the case of a system composed of N independent 
degrees of freedom, in the sense that the energy can be decomposed into a sum of N 
independent terms associated to each degree of freedom. In many cases, the density 
of states p{E) of such systems can be written in the form 

piE) « CiE - i^n.in)"^-! (18) 

for E close enough to i^inin (E > i?min), where C and rj are positive constants. For 
instance, in a perfect gas of N monoatomic particles, one has -Emin = and 77 = d/2, 
where d is the space dimension. 

We now determine the dynamical energy distribution P{E,t), again from the 
results of Sec. [31 The energy density p{E) satisfies Eq. with i' = rjN. Hence 
the dynamical energy distribution falls within the WeibuU class. The asymptotic 
distribution is then given by Eq. (IT6l) with a parameter 7 = — 1/j/ = —l/{riN), which 
is very small for large N. We know that the distribution ([T6| converges to the Gumbel 
law when 7 -^■ 0, thus for large N it becomes very close to the Gumbel distribution 

m- 

Hence, the connection to extreme value statistics presented in this paper allows 
one to determine in a straightforward way the dynamical energy distribution of the 
Barrat-Mezard model with an arbitrary density of states. In addition, it turns out 
that in many situations of interest, the dynamical energy distribution P(E,t) of a 
large system is actually very close to the Gumbel distribution, as can be seen on the 
above two examples. 
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5. Discussion: from aging to relaxation in dissipative models 

5.1. A toy model for undriven dissipative systems 

The Barrat-Mezard model has been defined within the framework of aging systems, 
in the zero-temperature limit. During the relaxation process, energy is progressively 
released to the zero-temperature reservoir. From a theoretical point of view, this 
situation is reminiscent of the dynamics of freely cooling dissipative systems, like 
undriven granular gases, or freely decaying turbulence. In such systems, the 
elementary degrees of freedom retained in the description are macroscopic entities 
(grains, etc.) or collective modes (for instance Fourier modes) coupled via dissipative 
interactions. In the example of the granular gas, energy is dissipated through collisions 
between grains, while in the case of turbulence, dissipation results from the viscous 
damping. 

In most realistic cases, dissipative events like collisions involve only a small 
number of degrees of freedom in the systems, and thus do not lead to large scale 
changes of the configuration of the system. However, in a mean-field spirit, one could 
model a dissipative event as a transition to a randomly chosen configuration with 
lower energy. Although very naive, this picture however leads to a reasonable scaling 
of the dissipated energy in each event, as a function of system size. Indeed, one 
expects in a gas composed of TV grains that the fraction of energy dissipated in a 
collision is proportional to since only a finite number of grains (typically two) 

are involved in the collision. Considering now the Barrat-Mezard model defined with 
the energy density p{E) oc E~^^, corresponding to a gas of particles, we compute 
the conditional average value of the energy E' reached after a single transition starting 
from energy E, yielding 

, ,^ \fdE'E'p{E') ^ , ^ 

{E')e = ^ = -TTT^ E. 19 

j^dE'p{E') ^ + 1 

It follows that the typical amount Ai^typ of energy dissipated to the zero-temperature 
heat bath reads, in the large N limit 

2E 

AEty^^E-(E'}E^^. (20) 

Hence, despite its mean-field nature, the Barrat-Mezard dynamics yields the correct 
scaling for the energy dissipated in a single event. This result suggests that some 
qualitative features of dissipative systems might be captured by the Barrat-Mezard 
model. 



5.2. Global fluctuations in granular gases 

Along this line of thought, and given the connection presented in Sec. [3] between 
extreme value statistics and the Barrat-Mezard model, it is natural to wonder 
whether the latter model could qualitatively account for the surprising emergence 
of distributions close to the Gumbel form in the relaxation regime of granular gases 
|14] . More precisely, it was found in numerical simulations of a two-dimensional freely 
cooling granular gas that the shape of the energy distribution P{E, t) is well-described 
by a Gumbel-like distribution [14 , also called BHP-distribution in this context [15]. 

However, the Gumbel distribution observed in the Barrat-Mezard model differs 
from that reported in granular gases |14] by a sign reversal. The former, being the 
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distribution of minimum values, has an exponential tail on the left side, while the latter 
exhibits an exponential tail on the right side. Hence, the simple cooling mechanism 
at play in the Barrat-Mezard model is not responsible for the appearance of the 
BHP-distribution in undriven dissipative systems. The reason for the appearance of 
Gumbel-like distributions rather seems to lie in the criticality of the dynamics due to 
the proximity of an instability threshold [111 [16] , as such distributions were originally 
reported in critical systems [15l [HI [18] . 

6. Conclusion 

In this short note, we have illustrated how the discrete-time Barrat-Mezard model 
at zero temperature, a simple model exhibiting 'entropic' aging, is quantitatively 
related to the extreme value statistics of independent and identically distributed 
random variables. Such a connection provides interesting insights into the long time 
statistics of the Barrat-Mezard model, showing that the dynamical energy distribution 
asymptotically corresponds to one of the limit distributions arising in extreme value 
statistics (Gumbel, Frechet or Weibull). This mapping has been established in the 
framework of the discrete-time model, but for times t ^ tq the time discretization 
becomes irrelevant, and the results for the long-time behaviour of the distribution 
P{E, t) should also hold for the continuous-time model. Note also that the connection 
presented here between glassy systems and extreme value statistics differs from that 
reported in [T2] . 

An interesting issue to be explored further would be the link to record breaking 
processes |19| . that have been put forward as an important mechanism at play 
in aging dynamics |20] . In the Barrat-Mezard model at zero temperature, all 
transitions occuring between configurations can be seen as record breakings, if one 
considers that random numbers are drawn at each time-step of the discrete-time 
dynamics, as proposed in Sec. |3l Another line for further investigation deals with 
possible connections between extreme value statistics and problems of sums of random 
variables. Such a connection has already been established in a different context 
[5TJ [221 [12] • There, the maximum in a set of independent and identically distributed 
random variables has been mapped onto a sum of non-identically distributed, and 
most often correlated, random variables. Such a mapping accounts for instance for 
the appearance of the Gumbel distribution in the statistics of the integrated spectrum 
in 1//- noise problems [UlUT]. In the Barrat-Mezard model, an analogous connection 
could possibly be made. Indeed, the energy reached at time t can be expressed as the 
sum of the energy at time t = and of all the energy jumps between times and 
t. However, the number n of terms in the sum (that is, the number of jumps) is not 
constant, and one expects n to be correlated in a non-trivial way to the total energy 
E, making the problem difficult to handle from the probabilistic side. Nevertheless, 
results obtained in a similar spirit for extreme values of sets of a non-constant number 
of random variables have been reported recently , suggesting that such statistical 
issues should be worth investigating. 
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